Emergence of rigidity at the structural glass transition: a first principle computation 
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We compute the shear modulus of structural glasses from a first principle approach based on 
the cloned liquid theory. We find that the intra-state shear-modulus, which corresponds to the 
plateau modulus measured in linear visco-elastic measurements, strongly depends on temperature 
and vanishes continuously when the temperature is increased beyond the glass temperature. 

PACS numbers: 61.43.Fs 



The shear-modulus is an unambiguous measure of the 
mechanical stability of materials. When one cools a glass- 
former below its glass transition, there appears a non- 
zero shear modulus on laboratory-accessible time scales. 
The understanding of the mechanism through which this 
rigidity emerges at the glass transition is a basic problem 
in condensed matter physics. 

A standard view on glasses is to regard them as 
very slow liquids with extremely high shear- viscosity [ij . 
Visco-elastic measurements show that super-cooled liq- 
uids and various soft-glassy materials behave as solids: 
The elastic modulus develops at low frequencies a 
plateau, which extends to lower and lower frequencies by 
lowering temperature or increasing density (see 0,01 a-nd 
references therein). Thus glasses aquire rigidity progres- 
sively. This feature is remarkably different from ordinary 
transitions from liquid to crystal where the rigidity ap- 
pears abruptly at the 1st order phase transition. 

Among various theoretical attempts, the so-called ran- 
dom first order theory (RFOT)II provides a useful work- 
ing ground to study the super-cooled liquids and glasses 
in a unified manner [5|. At the mean field level it is 
backed up by some microscopic approaches. On the one 
hand, it is intimately related to the mode-coupling theory 
(MCT) concerning the dynamics at relatively high tem- 
peratures 0. On the other hand, the so-called cloned liq- 
uid approach which combines the traditional liquid the- 
ory and the replica method allows one to compute ther- 
modynamic static quantities at lower temperatures 
This approach is currently the main first-principle ap- 
proach to studying properties of the glass phase. It has 
been so far limited to computing thermodynamic prop- 
erties, in particular the 'complexity' giving the entropy 
associated with the number of glass states. One of the 
major present challenges is to understand how the nu- 
cleation processes allowing to jump between glass states, 
which are not taken into account in the simplest RFOT 
scenario, can be included into this scheme. These pro- 
cesses are in particular crucial to explain why the mean 
field prediction of a dynamical transition at the dynami- 
cal (MCT) temperature Tc breaks down, and is replaced 
by a rapid increase of the relaxation time when the tem- 



perature gets close to Tc- Attempts in this direction in- 
clude the mosaic theory of Q and the study of long-range 
interactions [lo| . 

We shall extend the cloned liquid approach in order to 
compute the static shear-response of glasses. We identify 
the plateau modulus mentioned above by distinguishing 
intra-state and inter-state stress fluctuations. When ap- 
plied to mesoscopic samples, this approach predicts that 
the stress vs strain curve should have an intermittent 
behaviour. 

Models - We consider a system of N particles {i = 
1,2,..., N) at position = (x^, yi, zi) in the laboratory 
frame, which are interacting with each other via two body 



potentials, H{{vij}) 



and ' 



{Tij) , where 



|ry |. In order to study the rigidity against sim- 



ple shear deformation, we consider a system of particles 
in a container with two boundary walls which are normal 
to the 2-axis and separated from each other by distance 
Lz as shown in Fig[T] To impose a shear strain 7 on the 
system, we simply displace the top wall by an amount 
into i-direction. Then it is convenient to introduce 
a sheared frame with x',?/' and z' which are related to the 
laboratory frame as, (x, y, z) — {x' + z'7, y', z'). The vol- 
ume of the system V and the number density p = N/V 
remain constant under this shear. 
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FIG. 1: Schematic picture of the system under static shear, 
a) The mobile particles (open circles) are bounded by "ran- 
dom walls" (filled circle), b) Schematic mean-field picture of 
energy landscape and stress-strain curve 



Shear-modulus - The total free-energy i^(7) of the 
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system can be written as 
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with /3~^ = k^T being the inverse temperature. The 
ideal gas part of the free-energy is omitted because it 
does not change. Note that the integration over the in- 
terior volume of the container is taken using the sheared 
frame r' = {x',y',z') for which the integration range V 
is independent of the shear 7. 

Taking infinitesimal shear, the free-energy can be ex- 
panded formally as F(7) = F{0)+N{a)-f+^^lJ^+O{J^), 
where the shear-stress a and the shear-modulus ^ are: 
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where (...) denotes a thermal average evaluated with zero 
strain 7 = 0. We have introduced short hand notations 
like Xij = {xj — Xj)/rij, and the prime stands for dif- 
ferentiation. The 1st term (6) in the definition of fj, in 
Eq. ([2]) is called the Born term. It represents the instan- 
taneous response of the system against shear G(r = 0) 
(see below), which is finite even in liquids. The 2nd term 
is the correction term due to thermal fluctuations of the 
shear-stress. 

This type of fluctuation formula for the static elas- 
tic constants is well known [isl 25 1. In liquids, it is 



equivalent to the static limit of the Green-Kubo formula 
which relates the dynamic linear response against shear 
6a{t) = dt'G{t — t')j{t') to the shear-stress auto- 
correlation function {a{t)a{t')) by G(r) ~ /3((T(r)cr(0)). 
Linear visco-elastic measurements give access to the com- 
plex dynamical modulus G{uj) = ujG*{ijj). 

Boundary condition to shear - Let us pause here 
to discuss more explicitely the boundary condition. First 
of all, it is obvious that the boundary walls should not be 
strictly translationally invariant to exert shear on the sys- 
tem. On the other hand, we wish the system to maintain 
translational invariance at least on macroscopic scales. 
We thus assume that the walls are built from a quenched 
random configuration of particles, as shown in Fig. [1] so 
that the system keeps translational invariance in a sta- 
tistical sense. 

Because of the translational invariance at macro- 
scopic scales, the thermodynamic free-energy density 
limy^oo F{'^)/V is independent of7- This also happens 
in crystals, as discussed recently [ll[. However the def- 
inition of the shear modulus (and of the solid state) is 



through the linear response to a shear: the shear modu- 
lus in solids is non-zero because of the non-commutation 
of the small shear limit 57 —> and thermodynamic limit 
V ^ 00. The same phenomenon happens in glasses. This 
discussion has interesting consequences if one studies the 
deformation of mesoscopic samples on very small scales. 
We expect that the stress (tr) and shear-modulus ji of 
a single realization of the random walls will be non-zero 
even after the thermal averaging but fluctuate along the 
7-axis (See Fig. [T|). Only if one performs an average 
along the 7-axis will one recover the zero average value. 
Physically the breakdown of the commutation of the two 
limits means that elasticity theory fails. Thus elastic- 
ity and plasticity must emerge simultaneously in solids. 
We will argue that the plastic events can be viewed as 
changes of the relevant metastable states when one varies 
7. (see Fig.[T]). 

Shear on a cloned system - Let us now analyze 
the static response of glasses to shear. Taking the view 
of the RFOT, we suppose that there exist exponentially 
many metastable states a = 1,2,... with free-energies 
per particle /„. 

Our strategy is to consider a cloned system: m repli- 
cas (a — l,2,...,m) are forced to stay in the same 
metastable state, and we examine how the system re- 
sponds to a generalized shear such that each replica a is 
submitted to a different strain 7a. The total free-energy 
Fm of such a cloned system can be formally as 

m ^ ra 

F™(7) = ^^m(0)+A^^(a„)7a + - Ma67a7b+0(7') • 

a=l a.b=l 

(4) 

By clone symmetry the generalized shear-modulus [L^h 
can be written as 



(5) 



with 



W - ^/3EJ('^')a - ^o)l\P^ and 

Y.MlPc^-^Y.Mo.Pc.f\ where P„ 

is the thermal weight of the a-th 
metastable state at temperature T/m and (. . .)„ stands 
for a thermal average within the a-th metastable state, /i 
can be naturally interpreted as intra-state shear-modulus 
and jj, as the negative correction due to inter-state ther- 
mal fluctuations. It is easy to see that the physical shear- 
modulus ^ of a single system at temperature T/m can 
be obtained as /i = X^j^i Mob = ft + rnji. 

Physically the intra-state shear-modulus fi should be 
interpreted as the plateau modulus measured in linear 
visco-elastic measurements We expect it not to fluc- 
tuate between different metastable states. On the other 
hand, f, which is due to the inter-state fluctuations should 
be different on different realizations of the random walls. 
As we discussed before, the statistical translational sym- 
metry of the random walls requires the total modulus fi 
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to vanish on average; this imposes that, on average over 
the realizations of random walls, jl = —[i/m. 

Within the RFOT Q , the metastable states disappear 
at a dynamical transition temperature Tc(> Tk) pre- 
dicted by the MCT @. Thus the plateau modulus /t 
is positive only below Tc- 

Cage expansion of the shear-modulus - Within 
the cloned liquid theory (see details in Q) one assumes 
that particles in different replicas form molecules with 
a certain 'cage' size A, which plays the role of an order 
parameter that distinguishes the liquid phase A = oo 
from the glass phase A < oo. The system is considered 
as a liquid at an effective temperature T* = T /m* where 
TO* = m*{T) is determined for each temperature T by 
the stationnarity condition of the free energy. One finds 
that m* < 1 when T < Tk, where Tk is the Kauzmann 
temperature[l^, while to* sticks to the value 1 for larger 
temperatures. In many cases the behaviour of to* is well 
approximated by m* ~ T/Tk 0, S|- It is convenient to 
label the molecules as i = 1,2, ... ,N and to write the 
position of a particle as r^ = + uf where is the 
center of mass position of the molecule and uf describes 
the displacement of the particle in replica a within the 
molecule (X^a"? ^ "^^^ '^^S® ^i^*^ ^ = {Nm{m — 
Si Sa<b(("i'~"i)^) assumed to be small enough 
to allow a small cage expansion. Then the fluctuations 
within the cage is characterized by {ufu^)cagc = ^2(1 — 
TTL5ab)^5ij, where m is a component of u. If the cage size 
A changes discontinuously from a finite value to c» at 
Tc, as predicted by the MCT the cage expansion can 
work in principle right up to Tc from below [9| . 

Using the above prescription, we have computed the 
shear-modulus up to 1st order in the cage expansion 
The intra-state (or plateau) modulus is obtained as 



A = (fc>* - Jl + (J2 + ^3)(l - m) 



(6) 



where (...)* is a thermal average at temperature T* and. 
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i<j k<l 



jij diiu uij arc the sunimands in Eq. ([3]) and 



where aij and h 
stands for a connected correlation function at T* . In the 
derivation of the above result, we used the fact that the 
shear-modulus is zero in the liquid (&)* - l3*{a'^)c* = 0. 

The remarkable fact is that we can compute the 
plateau modulus at temperatures between Tr- and Tc. 
From Eq. ([6]) , one finds that in this regime (where to* = 
1), it takes the simple and suggestive form ft ~ (6) — Ji. 
The physical interpretation of this result is very simple. 



Instantaneous shear modulus — 
Intra-state shear modulus 



FIG. 2; Shear modulus of a binary soft-sphere system, com- 
puted from the cloned liquid approach. The bottom curve 
gives the plateau shear modulus, it is positive below Tc — 
0.22. 



On time scales shorter than the a-relaxation time, the 
stress-field is essentially frozen in time. There the only 
appreciable fluctuations are those associated with the (3- 
relaxation. The term Ji represents the strength of stress 
fluctuations due to these processes. 

A test case: a binary soft-sphere system - To 
test our scheme we performed an explicit computation 
of the shear modulus of the standard binary mixture of 
particles with soft-core interactions [16[ . The Kauzmann 
temperature of this system is Tk — 0.14 [§| while the dy- 
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namical (MCT) transition temperature is Tc ~ 0.22 
In Fig. [2] we show the result for the case of density p = 1. 
To evaluate to* and the radial distribution function 5*(r) 
at T* , we performed the cloned liquid computation using 
the binary HNC approximation [8| . 

The evaluation of various terms in Eq. (|6]) is done as 
follows. The Born term and J2 involve only 2-point func- 
tions so that they can be evaluated easily using g*{r). 
We evaluate Ji by Ji = ~2^ / d^rg^{r)P*\Va\'^{r) - 
2^/ d3rid3r2g(ri)g(r2)g(ri2)/3*V(T(ri) • Va{r2) where 
we made a chain approximation in order to approximate 
the three-point correlation function by a product of two- 
point terms, an approximation which is reasonable at 
high densities. The evaluation of J3 involves a connected 
4-point function which is expected to be smaller than the 
other terms and we have neglected it at present. 

As shown in Fig [2l the plateau-modulus /x strongly 
depends on the temperature. Remarkably, it continu- 
ously crosses at a temperature very close to Tc deter- 
mined by direct numerical simulations [l7j . Furthermore 
the cage size A is found to be still very small (10-2) at 
this temperature, which justifies our use of the first or- 
der small cage expansion. Our result of a plateau shear 
modulus emerging continuously below Tc disagrees with 
the conventional MCT , which predicts that the shear- 
modulus jumps discontinuously to a finite value at Tc 
from the liquid side. Our result means that the density 
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field is frozen at Tc as the MCT predicted, but the sys- 
tem is just marginally stable there [l^ . a picture which 
is consistent with the energy landscape picture of the 
RFOT [l8l - [2^ . For the visco-elastic measurements, this 
continuous transition suggests a power law behaviour 
G"(w),G"(w) oc Lo^. 

Intermittency of static shear response - Our re- 
sults have a natural interpretation within a mean-field 
picture. They suggest the following "intermittent" na- 
ture of static shear response below Tk at mesoscopic 
scales such that the system size N is large but finite. 
Within mean field, the parameter m = m*{T) ~ T/Tyl 
obtained in the cloned liquid approach is naturally in- 
terpreted as the Parisi parameter of the 1 step replica 
symmetry breaking (IRSB) Ansatz for the glass phase. 

As shown in Fig. [1] the interpretation at the mean 
field level is that of a free-energy landscape ^(7) which 
may be viewed as sequence of parabola with curvature /t 
(plateau-modulus) along the 7-axis, matching with each 
other at yield points [23[. Below Tk, the static response 
to shear is dominated by intra-state response with occa- 
sional inter-state response when passing the yield points. 

This picture is analogous to the mesoscopic response in 
mean- field spin-glass models [2J|. Here 7 plays the role 
of the external magnetic field h in spin glasses, which 
exhibit step-wise increase of magnetization m(h) along 
/i-axis. The drops of the stress passing the yield points 
corresponds to steps of the magnetization. At a given 7, 
each metastable state has a random free-energy and a 
random stress {<j)a so that the increase of 7 induces level 
crossings between low-lying states. 

The distribution of the stress may be modeled by a 
Gaussian distribution with zero average and variance 
A/^/N. From the correspondence with the spin-glass 
problem [24| we expect the typical spacing between the 
yield points to scale as Sjs 7k/(A-\/]V) and the 
width of thermal rounding of the yield points to scale 
as (57w ^ T/{A^/N). Here the parameter A is fixed 
as /3A^ = /t/rn in order to satisfy the condition that 
the total shear-modulus, including the inter-state shear- 
modulus, becomes zero on average. At low temperatures, 
if we choose a value of 7 randomly, most of the time we 
will observe the plateau modulus /t which is positive, and 
occasionally, with probablity Sjv//S% ^ T/Tk, we will 
find a negative shear-modulus. 

Discussion - Our computations predict a non-zero 
plateau-modulus /i at all temperatures below the dynam- 
ical transition ternperature Tc, including in the low tem- 
perature regime [25|. They also give a natural way to 
compute this dynamical transition temperature within 
the cloned liquid theory, offering an alternative to the 
MCT computation. This plateau modulus should be ob- 
servable dynamically on time scales smaller than the a 
relaxation time Tq,. Therefore one expects it to be seen, 
on all laboratory time scales, at all temperatures below 
the glass transition temperature (where the a relaxation 



time becomes larger than lO^s). 

The prediction of intermittent shear response in meso- 
scopic samples should also be amenable to experimental 
tests. It is supposed to take place even at temperatures 
higher than Tk at the length and time scales of the so- 
called mosaic states proposed by the RFOT [3, Q because 
each mosaic is subjected to a random pinning field pro- 
vided by surounding mosaics. 
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